Electronic Raman scattering in superconductors on triangle lattices 
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In superconductors on triangle lattices, we show that the electronic Raman spectra in the non- 
resonant channel are independent of the light polarization configuration in the scattering. Therefore 
the light polarization configuration is no longer a judicious choice to probe the anisotropic nature 
of the pairing wave function (as in the case of square lattices). However, we show that the detailed 
energy-shift dependence of the spectrum and the efi'ects of impurities provide useful information 
that can help identify the possible pairing symmetry experimentally. These results are relevant to 
the newly discovered Naa;Co02 • 2/H2 superconductors. 



The electronic Raman scattering in the so-called non- 
resonant channel (without invoking inter-band scatter- 
ing) measures the fluctuations of the effective density 

operator = Y.-ka^^^t+^a^^'^ of q — > 

(In the scattering event, q is the momentum transfer 
from the lights, which is vanishingly small compared to 
the Fermi momentum.) Here 7k is the Raman vertex 
that describes the diamagnetic second-order coupling be- 
tween electrons and lights, 7k = (e^ • Vk)(es • Vk)ek, 
where and are the polarization directions of the in- 
cident and scattered lights, and Ck is the normal state 
dispersion. The Raman response function w), with 
Lj — LJi ^ LOg being the energy transfer in the scattering, 
is the Fourier transform of the retarded Green's func- 
tion x(q, i) = — i([/Oq(t), p_q(0)])^^(^). The Raman in- 
tensity (i?/) is eventually given by — Im[x(q 0,w)]. 
Theoretically, x(Qi ^) can be conveniently obtained from 
the analytical continuation of the Matsubara propagator 
X(q, iuJn ^ + iO+). 

The electronic Raman response vanishes in ideal met- 
als because of lack of phase space for finite-energy and 
zero-momentum particle-hole excitations, but is finite in 
a superconductor because a Cooper pair can be bro- 
ken with no change in the center-of-mass momentum. 
By adjusting the incident and scattered photon polar- 
izations, through the Raman vertex function Raman 
scattering may projects out various parts of the quasi- 
particle excitations in the momentum space, and is there- 
fore highly useful to probe the angular dependence of 
the possible anisotropic pairing wave function on the 
Fermi surface. This has been applied very successfully 
in cuprate superconductors. [1] Because of the square lat- 
tice of the copper oxide plane, there are two Raman- 
active anisotropic channels, i.e., the Big and i?2g chan- 
nels. These channels can be conveniently picked up by 
adjusting the polarizations of the incident and scattered 
lights. The fact that the Raman responses are differ- 
ent in these channels is a proof that the pairing wave 
function is anisotropic, and is consistent with the d-wave 
symmetry. Furthermore, the detailed energy-dependence 
of the response yield useful information on the dynamics 
of the quasi-particles. It is therefore natural to explore 
how electronic Raman scattering can be used in super- 



conductors on other lattices. 

In this paper, we shall concentrate on triangle lat- 
tices, motivated by the discovery of the Naa;Co02 • 2/H2O 
superconductor, [2] which does have triangular C0O2 lay- 
ers, separated by thick insulating layers of Na+ ions and 
H2O molecules. The two-dimensional layered structure 
is similar to the copper oxide and it is conjectured that 
a similar superconductivity mechanism may be at work. 
In particular, it is argued that this is a doped Mott in- 
sulator that would be described by the resonant- valence- 
bond (RVB) theory. [3] In such a scenario, the pairing 
symmetry would be d^2_y2 -j- idxy with a full energy gap 
for quasi-particle excitations. On the other hand, given 
the experimental fact that the parent compounds exhibit 
intra-plane ferromagnetic correlations, px + ipy- [4] or /- 
wave [5] pairing symmetries are also proposed as possible 
candidates. For these triplet channels, the p-wave pair- 
ing has a full gap, while the /-wave pairing have nodal 
lines in the gap function. Finally, if the pairing is due 
to lattice vibrations, an s-wave gap function is to be ex- 
pected. As an important step toward understanding the 
pairing mechanism, one is therefore asked to discriminate 
experimentally one from the others in this list of possible 
pairing symmetries. 

Existing experimental data are still controversial. 
From nuclear-magnetic-resonance (NMR) measurements, 
Kato et al proposed p-wave triplet pairing, [6] while 
Kobayashi et al conjectured s-wave singlet pairing. [7] 
Nonetheless, from the coherence peak in the temperature 
dependence of the spin-lattice relaxation rate 1/Ti both 
groups agreed that the gap has no nodal lines. This is to 
be contrasted to the evidence of nodal lines in nuclear- 
quadruple-resonance (NQR) data [8] and in the specific 
heat measurement. [9] The main difficulty for all measure- 
ments is the control of the sample quality, and it is quite 
likely the reason behind the inconsistency of the experi- 
mental results. While this asks for better sample quality 
to resolve the inconsistency, in this paper we propose 
electronic Raman scattering as additional but important 
probe to the pairing symmetry. We note, however, that 
the low transition temperature (about 4.5K) may impose 
a difficulty for the measurement of the pairing gap due 
to the laser heating effect. Our hope is that the technical 
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difficulty could be removed and that samples with higher 
transition temperatures could be achieved. 

The structure of the rest of the paper is as follows. 
We first analyze the Raman response in the clean sys- 
tem, reaching the result that the response is independent 
of the polarization configuration of the scattering. We 
show the result by explicit calculation of a lattice model, 
and substantiate by a proof in the continuum limit. The 
Raman response in the clean system can only tell whether 
there is a full gap on the Fermi surface. Second, the effect 
of impurity is taken into account. We find that the im- 
purity effects do not induce Raman absorption below the 
pair breaking peak (PBP) for s-wave pairing, while they 
do so for px + ipy- and — y-^ + idxy-w&ve pairing. We 
also find a universal Raman response in the case of nodal 
/-wave pairing in the limit of small Raman energy shift 
for small scattering rates. Combined with a probe to the 
singlet versus triplet pairing, these features enables one 
to identify the pairing symmetry unambiguously. 

We start from the effective mean field Hamiltonian, 

H = Y1 ^kC^Ck. + E(^k4^Clki + h.c), (1) 

k(7 k 

where ek = — 2t^^^j^cosk ■ 5i — ji. Here 61,2,3 = 
(1,0), (1/2, V3/2) and (-1/2,^3/2) are three unit vec- 
tors on the triangle lattice, /i is the chemical potential. 
For the Na2;Co02 -2/1120 superconductor, band structure 
calculation reveals a large hole pocket around the T point 
in the momentum space, [10] it is therefore consistent to 
take t < in the dispersion. We also tune the chem- 
ical potential so that there are 1.35 electrons per site. 
Ak is the pairing gap function. For the s-, Px + ipy, 
(ij.2_j,2 + idxy- and /-wave pairing symmetries, the gap 
functions we will consider are, respectively. 



limit of q = can be written as, 

']y X] Tr[7k(T3G(k,ia;„)7k(T3G(k,iw„ -|-if„)],(6) 



k,iai„ 



where A'' is the number of lattice sites, and G is the mean 
field Green's function given by 

G~^(k,ia;„) = iw„(To - eucrs - (Ak(T+ -|- h.c). (7) 

Here a;„ and Vn are Fermionic and Bosonic Matsubara 

frequencies, respectively, (Tq is the 2x2 unitary matrix, 
ci,2,3 are the Pauli matrices, and = {ai+ia2)/2. The 
summation over Matsubara frequencies can be completed 
exactly. After the analytical continuation, we get the 
bare retarded Raman response function. 
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(8) 



where £^k = v^ej^ + |AkP. In the limit where 7k ~ 1 the 
response function is equivalently a bare charge density 
susceptibility in the long wave length limit. However, in a 
charged system, the spectral weight of this susceptibility 
is pushed up to the plasma frequency. While this fact can 
be taken into account by a random-phase-approximation- 
like approach, the effect of Coulomb screening on the 
Raman resposne can be taken into acc;ount by simply 
projecting away the symmetrical part as follows[ll] 
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Ak = A^cosk- Jj; 

i=l 

3 

A^ = A ^ sin(k • Si) exp(i6»i); 

i=l 
3 

Af = A ^ cos(k • 6i) exp(2i6»i); 

i=l 
3 

Aj^ = A ^ sin(k • Si) exp(3i^i). 



i=l 



(2) 
(3) 
(4) 
(5) 



where ^1,2.3 = 0, 7r/3 and 2tt/3 are the azimuthal angles 
of (5i,2,3, respectively. In the following calculation we take 
the gap amplitude as A = 0.1|t| for illustration purpose. 
The conclusion does not depend on this parameteriza- 
tion. We note that the gap on the Fermi surface is finite 
for the cases of s-, Px+ipy, and dx2-y2 +idxy-'wa.ve cases, 
and it has nodal points for the /-wave pairing. The bare 
Matsubara propagator for the effective density pq in the 



which is applicable below the plasma frequency. Here we 
suppressed the frequency arguments for brevity, and x°i 
and Xii are obtained by replacing 7k by 1 for one or two 
of the Raman vertices in x°^. 

To proceed, we write according to the given definition 
the three typical Raman vertex functions in the xx, yy 
and xy polarization configurations as follows. 



XX r>j./ 1,1 ^'X ^/Sky 

-ll^ = 2t{coskx +2^°^Y 



7r 



3t cos — cos 
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7k 



XV ATj. • ^X . ^/Sky 

w = -v3tsm — sm — — 



(10) 

(11) 

(12) 



The vertex function for an arbitrary polarization config- 
uration can be expressed in terms of the above as. 



7k = cos ( 
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FIG. 1: The electronic Raman spectra for px + ipy-wscvc pair- 
ing in the (a) xx, (b) yy and (c) xy polarization configura- 
tions. The dashed (solid) lines denote the bare (screened) 
spectra, (d) is the screened spectra in the three configura- 
tions plot together. Note that in (c) the dashed line falls on 
top of the solid line and is therefore invisible. 



where (j)i_s are the aziniuthal angles of the polarizations 
of the incident and scattered lights, respectively. 

With the above ingredients, we are now able to cal- 
culate the Raman spectra (the imaginary parts of Xj-y 
and X77) explicitly. We take the zero temperature limit 
in what follows. Fig.l is the spectra for the px + ipy- 
wave pairing in the xx (a) , yy (b) and xy (c) polarization 
configurations. The dashed (solid) lines denote the bare 
(screened) spectra. The PBP occurs near uj = 4A. This 
is because the effective gap on the Fermi surface is near 
2A, and the major contribution to the summation over 
k in x° is from the region near the Fermi surface. It is 
clear that while the bare spectra depends on the polar- 
ization configurations, the screened ones do not, which 
is made clearer by plotting them together in Fig. 1(d). 
This is true for any other arbitrary configurations (not 
shown here). Experimentally it is the screened spectra 
that are measured. We therefore conclude that the Ra- 
man spectra on triangle lattices are independent of the 
polarization configuration. 

In order to understand the above result, it is useful 
to resort to symmetry analysis. In general, the Raman 
vertex function can be decomposed in terms of the basis 
function of the irreducible representations of the point 
group. For the system under concern, there are only 
two Raman-active irreducible representations, i.e., the 
one-dimensional Aig and two-dime nsional Eg represen- 
tations. The corresponding basis function for Aig is 



7,0 \/2>ky 

7q = coskx -I- 2 cos Y cos ^ , 



(14) 



and that for the Eg representation is 7e = (7ei,7e2) with 



/t: • kx . ^/Sky 
el = V3sm ysm^-^; 



kx y/Sky 

7e2 = cos kx - COS — COS ^ " . 



(16) 



Symmetry arguments guarantee that the Raman re- 
sponse is also decomposed into these two channels. Fur- 
thermore, since 7^ has the full symmetry of the system it- 
self, the Raman response in this channel will be screened. 
Therefore, Eg is the only channel for electronic Raman 
response, and this implies that at least the line shape of 
the Raman response should be independent of the polar- 
ization configuration as illustrated in Fig.l. It remains to 
be shown that even the amplitude of the Raman spectra 
is independent of the polarization configuration. We first 
verify the claim by taking advantage of the typical result 
in Fig. 1 , and then present a more rigorous proof by going 
to the continuum limit. 

In terms of the basis functions, the three typical Ra- 
man vertex functions can be rewritten as. 



7'''' =*(7a+7e2); 
7^^=t(7a-7e2); 

7"^ = -t7el- 



(17) 
(18) 
(19) 



Since the 7q component does not contribute to X77 due to 
screening, Fig.l proves that the 7ei,e2 components con- 



tribute equally to the X771 or X7ei' 



X7e; 



Now 



consider an arbitrary vertex which can be rewritten as 

7 = t(c0S - sin 07e2), (20) 

where cf) = <j)i + 4>f , and we dropped the 7a component, 
anticipating the screening effect. Consequently, 



X77 = COS^ <t>Xjellel + sin 0X7e27e2 = Xle 



(21) 



where a contribution from the crossing term 7ei7e2 van- 
ishes by symmetry. 

To close the argument, we now prove X7ei7ci ^ Xieii^i 
for all of the pairing symmetries listed previously in 
the continuum limit. Since the screening effect is ab- 
sent in the anisotropic Eg channel, it suffices to prove 
x£ ,™ , = x!^ ,-7 Near the Fermi surface, one can write 
approximately but without loss of symmetry (up to triv- 
ial common factors). 



7ei ^ sin 29] 7e2 ^ cos 26, 



(22) 



where 9 is the azimuthal angle of the Fermi momentum 
kjr approximately on a circle. On the other hand, the 
gap functions can be approximated as 



~ A^, 

Ai? exp(i6'), 
A^ ~ Af exp(2i6i), 
A^ ~ AFsin(36»), 



A? 



(23) 

(24) 
(25) 
(26) 



where Ap is the effective gap amplitude on the Fermi 



(15) 

surface for the specific cases. It is seen that jAg 
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either a constant (for s-, px + ipy- and dx2^y2 + id^y- 
symmetries) or six-fold symmetric (for /-wave symme- 
try) on the Fermi surface, whereas ^2 = (l±cos40)/2 
which is four-fold symmetric. Plugging these into the ex- 
pression for x° and going over to the continuum limit, 
one concludes that by symmetry only the constant part 
in 7g]^ g2 contributes to the response function. This com- 
plete the prove. It is interesting to point out that in 
the case of da;2_j^2-wave pairing on a square lattice, jei,e2 
becomes the basis functions for the B2g and Big rep- 
resentations. In this case, [A^ip is four-fold symmetric, 
which can resonate with the corresponding component in 
7p]^ p2 ; and this is why in that case the Raman spectra is 
different in Big and i?2g channels. 

Having convinced ourselves that the Raman spectra on 
a triangle lattice do not depend on polarization config- 
uration for all the relevant pairing symmetries, we now 
look into the energy dependence of the spectra. Fig. 2 
shows the screened Raman spectra for the four pairing 
symmetries under concern. As anticipated, the spectra 
have a threshold for the fully-gapped (a) s-, (b) p^+ipy-, 
and (c) (ia;2_j,2 + id xy-wavc pairing, whereas it is linear at 
low energy for the nodal (d) /-wave pairing. The PBP 
in Fig. 2 occur at different energies in units of A. This is 
again due to the fact that the effective gap amplitude Ap 
on the Fermi surface differs in the cases considered, even 
though the same A is used in the gap function. We also 
note that the line shape for the d-wave case is slightly dif- 
ferent from the s- and p-wave cases. This is because the 
d-wave gap value under concern is close to the van- Hove 
singularity of the normal state density of states (DOS), 
which contributes a secondary peak near the threshold. 
We do not expect it to be a generic feature. In summary, 
in the clean system considered so far, the Raman spec- 
tra can only tell whether the pairing gap is nodal or not. 
However, impurities are always present in real systems, 
and it is pending to see whether impurity effects can help 
differentiate the three fully-gapped pairing cases. This is 
to be anticipated as impurities may induce low energy 
(Andreev) states for nontrivial pairing, even if it is fully 
gapped in the clean system. This will be discussed in the 
following. 

The effect of impurities can be conveniently accounted 
for by the coherent potential approximation (CPA). The 
modified Green's function is determined by 

G~^(k, iojn) = iconcro - ekO-3 - (Akcr+ + h.c.) - T,{oJr(?7) 

where E(a;„) = So(a;„)c7o + Si(ci;„)(7i. Henceforth we 
apply the continuum limit and assume particle-hole sym- 
metry near the Fermi surface so that the (73 component 
of S vanishes. We take A£. to be real for definiteness. On 
the other hand, we anticipate that the ai^2 (or a^) com- 
ponents of S vanishes for pairing symmetries other than 
the s-wave, as to be verified hereafter. The self-energy S 
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FIG. 2: The screened electronic Raman spectra for (a) s-, (b) 
Px + iPy-, (c) d^2_y2 + idxy-, and (d) /-wave pairing. 



is determined by the self-consistent condition 

r 



9 - cuzg ^0-3 ' 



(28) 



where F = rii/nJVo is the scattering rate with and 
Mo being the density of impurities and the DOS at the 
Fermi energy, respectively, c = cot 5 with S being the 
scattering phase shift, and g is the normalized on-site 
Green's function 



irAfoJ (27r)- 



G(k, iwn). 



(29) 



It is straightforward to verify that the cts component of 
g is absent by particle-hole symmetry, and so are the 
CTi_2 components in the case of non-s-wave pairing, jus- 
tifying the a priori assumptions. Similar self-consistent 
equations holds for the retarded and advanced Green's 
functions, and are used in the realistic calculations via 
the spectral representation, 



1 /■ G'^(k,w)-G«(k,c 



ILJn — 



-dw. (30) 



This is substituted into Eqs.(6) and (9) to calculate the 
Raman spectra in the presence of impurities. [12] Vertex 
corrections have been ignored, and are not expected to 
affect the conclusions significantly. 

We first discuss the effect of unitary hmit (c = 0). 
Fig. 3 (a) and (c) show the Raman spectra for s- and 
Px + ipy-wavc pairing under various values of the scatter- 
ing rate F. For the s-wave pairing, the Raman spectra 
remains forbidden below the PBP at a; = 2 Air. Impu- 
rity scattering merely changes slightly the spectra above 
the peak. This is easily understood from the behavior 
of the DOS plot in Fig.3(b), where impurity does not 
induce intra-gap quasi-particles, consistent with the An- 
derson theorem. The slight change in the Raman spectra 
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FIG. 3: The electronic Raman spectra for (a) s- and (c) Px + 
ipy-wave pairing in the unitary hmit for various values of the 
scattering rate. The corresponding DOS are shown for (b) s- 
and (d) px + ipy-wave pairing. 



is due to the slight change in the DOS above the gap en- 
ergy. The situation is quite different in Fig.3(c), where 
the Raman spectra is induced below the PBP. This can 
again be understood from the DOS induced by impuri- 
ties in Fig.3(d). While there is a full gap in the clean 
limit, an impurity band is induced near the Fermi energy 
for small T. This is due to the nontrivial intrinsic phase 
of the Px + ipy-wa.ve pairing function. During impurity 
scattering, the quasi-particles experience a change of the 
intrinsic phase. It is shown elsewhere that by such an 
effect a single unitary impurity can induce an Andreev 
bound state, [13]and the bound states from a distribu- 
tion of impurities can therefore form an impurity band. 
The scattering of quasi-particles by the light from the 
impurity band (below the Fermi energy) to the gap edge 
(above the Fermi energy) contributes most significantly 
to the Raman absorption below the PBP. This also ex- 
plains the hump with a pseudo-threshold at uj = Ap. For 
very large F, the intra-gap is filled with induced DOS and 
accordingly the Raman spectra varies smoothly below 
2Af. Wc add that the Raman spectra for d^2_y2 +idxy- 
wave pairing behaves similarly to the px + ipy-wave pair- 
ing under a similar mechanism, and are not shown here 
for brevity. 

Fig. 4(a) shows the Raman spectra in the unitary limit 
for /-wave pairing. Although the PBP is lowered with 
increasing scattering rate, the low energy part is hardly 
changed. This implies a universal low energy Raman 
spectra in this case, to which we shall return shortly. The 
DOS is plot in Fig. 4(b). In the absence of impurities, 
the DOS is linear at low energy, and impurities again 
induce low energy states. However, the impurity band is 
different from that in the px + ipy-wavc case, in that this 
band connects smoothly to the higher energy DOS. This 
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FIG. 4: (a) The electronic Raman spectra for /-wave pairing 
in the unitary limit, (b) The corresponding DOS. 



explains the absence of new humps below the PBP in the 
Raman spectra. The universal low energy Raman spectra 
deserves further scrutiny. For a low scattering rate F, 
this can be proved exactly. Prom the CPA self-consistent 
equation the retarded self-energy T,^'^ = Hrjao in the 
low energy limit, where 77 ~ y/FAp. If this is substituted 
into Eqs.(30) and (6), we find after some algebra that for 
LU <^ Ap the Raman spectra is given by — Imx(w„ U! + 
«0+) = auij"^ /vpVA, where a is a constant, 7^ = {'y^)p 
is the average square of the Raman vertex on the Fermi 
surface, vp is the Fermi velocity and = |VkAk|k=k„ is 
the gap velocity at the six nodal points k„ on the Fermi 
surface. The result is tied to the existence of nodal quasi- 
partic;les. The universal behavior in the Raman spectra 
is reminiscent to the universal conductivity discovered 
in the literature for d^^s-yS-wave pairing in the cuprate 
superconductors. [14] 

We have also considered general values of the inverse 
scattering strength c. The conchision does not change for 
the s-wave case. The Raman spectra for px + ipy-wave 
case is shown in Fig. 5(a) for various values of c. Fig. 5(b) 
is the corresponding DOS. It is seen that near the Born 
limit, the center of the impurity band shifts away from 
the Fermi energy, causing a corresponding change in the 
Raman spectra. The Raman spectra for /-wave case is 
shown in Fig. 5(c), and the DOS in Fig. 5(d). Here no 
impurity band is formed in the Born limit, and the zero- 
energy DOS increases gradually with decreasing c. Inter- 
estingly this does not cause a change in the low energy 
limit of the Raman spectra. The reason is that for any 
values of c, as long as the resulting self energy is indepen- 
dent of energy in the low energy limit, the mechanism of 
universal Raman spectra applies as described above. 

In summary, we developed a theory of electronic Ra- 
man scattering in superconductors on triangle lattices, 
with possible application to the cobalt oxide supercon- 
ductors. The Raman response is shown to be indepen- 
dent of the light polarization configuration. However, in 
the case of s-wave pairing the absorption is forbidden 
below the PBP in the clean system as well as in the pres- 



6 



- - c=0 

1 

100 


(a) 


D 2 4 


- - c=0 

0.2 J 

100/ 

J/ 


(c) 



c " 

=3 

€ 4 

_co 

w 2 
O 

Q 



(b) 



-3 -2 



-1 
co/A^ 



1 2 3 




2 

co/A^ 



FIG. 5: The electronic Raman spectra for (a) p^. + iPy- and 
(c) /-wave pairing with F/Af = 0.1 and various values of c. 
The corresponding DOS are shown for (b) "Px + vpy- and (d) 
/-wave pairing. 



ence of impurity scattering; in the cases of 'Px + ipy and 
dx^-y^ + idxy-'weLVC pairing the absorption below the PBP 
is forbidden in the clean system, but can be induced by 
impurities; and finally in the case of nodal /-wave pair- 
ing there exists a universal low energy Raman response 
in the low energy limit for weak scattering rates. These 
features, combined with another probe of singlet versus 
triplet pairing, can be used to identify the relevant pair- 
ing symmetries unambiguously. 
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